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PREFACE 


GENERAL:  This  publication  is  one  of  a  series  of  Reference 
Publications  on  achi events  related  to  the  fields  of  mapping, 
charting,  and  geodesy,  and  their  related  arts  and  sciences. 

Each  Reference  Publication  is  written  by  a  Defense  Mapping 
Agency  Aerospace  Center  technician  qualified  by  training  and 
experience  to  contribute  knowledge  and  technology  to  the 
selected  subject. 

PURPOSE:  To  contribute  techni^l  information  to  the  field  of 
geodesy  by  describing  the*  computational  methods  used  ty  DMAAC 
to  determine  1°  x  1°  moan  free-air  gravity  anomalies  und  their 
accuracies. 

DISCLAIMER:  The  fi.'^'  .js  cf  this  report  are  not  to  be  construed 
as  the  official  Defei.ae  Mapping  Agency  position  unless  so  designated 
by  other  authorized  documents. 

REPRODUCTION:  This  publication  does  not  contain  copyrighted 
material,  nor  is  a  copyright  pending.  Reproduction  in  whole  or 
in  part  is  permitted  for  any  purpose  of  the  United  States 
Government. 

DISTRIBUTION:  General  public  release  by  the  National  Technical 
Information  Service,  U.S.  Department  of  Commerce,  Springfield, 
Virginia  22151,  is  approved.  Qualified  requesters  may  obtain 
copies  from  the  Defense  Docunentation  Center. 
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ABSTRACT 


This  publication  contains  the  explanations  of  the  methods  used 
1n/?!!P“tin9  the  !°  x  ]°  mean  gravity  anomaly  values 

™d  73*0Q02CCiraCieS  t0  be  published  Reference  Publication 

.  The  information  contained  in  this  report  is  in  two  parts.  Part  I 
identifies  each  computational  method,  gives  a  written  definition,  and 
funnshes  the  equations  with  a  sample  computation.  When  end  how  the 

aro^aH/)Sh0U  j  benus!d  and  any  strengths  and  shortcomings  of  the  method 
are  addressed.  Part  II  identifies  the  method  used  to  arrive  at  an 
accuracy  value  for  each  computational  method,  gives  a  written  definition 
with  associated  equations,  and  provides  a  sample  computation.  Short¬ 
comings  of  the  methods  are  stated. 
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COMPUTATIONAL  METHODS  '-OR  DETERMINING  1°  x  1°  MEAN 
GRAVITY  ANOMALIES  AND  THErR  ACCURACIES 

INTRODUCTION 

In  the  ideal  situation,  a  mean  free-air  gravity  anomaly  CP*f) 
for  an  area  would  be  obtained  from  an  infinite  number  of  well 
distributed  gravity  observations.  There  are  no  areas  where  the 
density  and  distribution  of  gravity  observations  permit  this  ideal 
situation.  In  most  cases,  the  gravity  observations  are  either 
concentrated  in  local  areas  or  along  transportation  routes  so  that 
it  is  necessary  to  compute  the  mean  free-air  gravity  anomaly  for 
an  areo  from  gravity  observational  data  with  irregular  distribution. 

This  necessitates  using  various  methods  to  compute  A'jf. 

The  Ag^  values  are  computed  for  all  areas  of  the  world  where 
sufficient  gravity  data  is  available.  In  areas  where  there  is  a 
requirement  for  a  Ag^  and  there  is  sparse  or  no  gravity  coverage, 
the  Agf  is  determined  using  gravity-geophysical  correction  techniques. 
The  Ag^  values  represent  the  mean  values  of  the  anomalous  gravity  within 
the  surface  areas  of  various  required  sizes,  and  are  developed  directly 
from  observed  gravity  data  whenever  possible.  In  ocean  areas  the  Ag^.  is 
computed  from  point  free-air  gravity  anomalies.  Over  land  areas  it  is 
computed  using  point  Bouguer  gravity  anomalies  and  mean  elevations. 

All  gravity  data  held  by  the  DOD  Gravity  Library  is  referenced  to 
the  International  Gravity  Standardization  Net,  1971  ( IGSN  71),  and 
gravity  anomalies  are  computed  using  the  Gravity  Formula  for  the 
Geodetic  Reference  System  1967  (GFGRS  67). 
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All  computed  ifg^  values  and  supporting  Information  are  recorded 
In  punch  card  form  and  stored  on  magnetic  tape.  The  following  Is  a 
list  of  codes  which  Indicate  the  method  used  to  compute/predict  the 
1°  x  1°  mean  free-air  gravity  anomaly: 

3  -  Bouguer  anomaly  map  estimates. 

4  -  Free-air  anomaly  map  estimates. 

6  -  Geophysical  correlation. 

8  -  Interpolated  or  extrapolated. 

9  -  Mean  anomalies  received  from  other  organizations. 

A  -  Average  of  smaller  size  quads. 

M  -  Modified  average  free-air. 

MEAN  GRAVITY  ANOMALY  COMPUTATION 

1 .  Bouguer  Gravity  Anomaly  Map  Estimates 

In  this  method,  a  rectangular  grid  with  compartments  of  desired 
;  'ze  Is  overlayed  on  a  contoured  Bouguer  gravity  anomaly  map.  The 
mean  Bouguer  gravity  anomaly  GTgg)  Is  estimated  fc^  each  rectangular 
compartment. 

The  a gf  Is  computed  using  the  following  equation: 

l 

_  _ 

=  AgB  +  0.1119  Fm 

where: 

^  =  the  estimated  mean  elevation  of  the  '•ompartment  in  meters. 
0.1119  mgal/m  =  the  effect  of  the  Bouguer  plate. 

I 
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This  method  Is  used  In  areas  with  Insufficient  point  gravity  data 
In  the  DOD  Gravity  Library  for  direct  computation  provided  a 
Bouguer  gravity  anomaly  map  Is  available 
Example: 


Figure  1.  Contours  of  Bouguer  Gravity  Anomalies 


The  mj  of  the  Inset  block  is  10  mgals,  and  the  .nean  elevation 
B 

estimated  from  a  topographic  map  is  50  meters.  The  itg  is; 

£g  =  10  mgals  +  (0.1119  mgal/m)(50m) 

Ag^  =  10  mgals  +  5.6  mgals 
£g^  =  15.6  mgals. 
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The  Bouguer  gravity  anomaly  map  method  is  slow  and  prone  to  human 
error.  The  estimates  made  by  the  analyst  are  subjective,  and  since 
this  is  basically  a  manual  method,  blunders  are  possible.  Grids 
must  be  made  for  each  map  used  and  the  grid  size  varies  for  each 
latitude  band. 

2.  Free-Air  Gravity  Anomaly  Map  Estimates 

Ihis  method  is  mechanically  the  same  as  the  Bouguer  Gravity 

Anomaly  Map  Estimates  method.  A  rectangular  grid  with  compartments 

of  desired  size  is  overlayed  on  a  contoured  free-air  gravity  anomaly 

map.  The  Agf  is  estimated  for  each  rectangular  compartment.  This 

method  is  used  for  areas  with  little  or  no  point  gravity  data  in 

the  DOD  gravity  Library  if  a  free-air  gravity  anomaly  map  is  avail abl 
The  Agf  of  the  inset  block  in  Figure  2  is  20  mgals. 


Figure  2.  Contours  of  Free-Air  Gravity  Anomalies 
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This  method  has  all  of  the  disadvantages  of  the  Bouguer  Gravity 
Anomaly  Map  method  plus  the  disadvantage  of  steep  gradients  inherent 
in  the  dependence  of  the  free-air  gravity  anomaly  on  elevation. 
Special  attention  has  to  be  given  to  bathymetric  features,  such  as 
trenches,  island  arcs,  etc.,  in  ocean  areas  which  usually  have 
strong  fluctuating  gravity  anomaly  fields  that  may  not  be  reflected 
in  a  free-air  gravity  anomaly  map  because  of  the  poor  distribution 
of  point  observations.  The  use  of  this  method  is  restricted  to 
water  areas  and  flat  land  areas. 

3.  Geophysical  Correlation 

Geophysical  correlations  are  used  to  compute  mean  free-air 
gravity  anomalies  in  gravimetrically  deficient  land  and  water  areas. 
The  many  different  techniques  used  are  being  documented  in  "Gravity 
Correlation  Study:  Manual  of  1°  x  1°  Mean  Gravity  Anomaly 
Prediction  for  Continental  Areas,"  in  preparation.  These  values 
are  recomputed  and  replaced  using  another  technique  as  soon  as 
sufficient  gravity  data  becomes  available. 

4.  Interpolation  and  Extrapolation 

This  method  is  primarily  used  in  ocean  areas  with  insufficient 
point  gravity  anomaly  data  to  compute  a  Ag^.  The  point  gravity 
'  anomaly  can  generally  be  correlated  with  the  bathymetric 
features  in  surveyed  areas,  and  the  same  correlation  applied 
in  unsurveyed  areas  with  similar  bathymetric  features.  The 
systematic  relation  to  water  depth  is  generally  linear  and  the 
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different  bathymetric  features  (islands,  trenches,  sea  mounts,  etc.) 
have  a  characteristic  influence  on  the  point  gravity  anomaly. 

The  analyst's  knowledge  of  the  expected  behavior  of  the  gravity 
field  is  used  to  aid  in  contouring  or  to  continue  contours  with 
respect  to  bathymetric  features.  The  Ag^  is  then  estimated  from 
the  contours. 

Example: 


Figure  5.  ExtrapolateJ  Anomaly  Field 
The  diagrams  illustrate  the  general  method  of  extrapolating  a  gravity 
anomaly  field  Starting  with  point  values  from  observations  (Figure  3) 
and  the  bathymetric  information  (Figure  4),  the  gravity  anomaly  field 


7 


for  the  quad  is  extrapolated  (Figure  5)  using  a  linear  relationship 
between  the  point  gravity  anomaly  and  bathymetry,. 

This  method  cin  be  used  on  land  with  Agg  as  point  data  instead 
of  the  Agf  because  of  the  smoother  field.  The  method  is  at  best 
an  approximation,  and  is  basically  a  manual  process.  These  values 
are  re'jmputed  and  replaced  using  another  technique  as  soon  as  sufficient 
gravity  data  becomes  available. 

5.  Mean  Anomalies  Received  from  Other  Organizations 

These  values  may  be  mean  Souguer  or  free-air  gravity  anomalies. 

The  methods  used  for  competing  the  mean  gravity  anomaly  (Ag)  values 
and  their  accuracies  vary.  Many  times  the  methods  used  are  either 
not  documented  or  not  available.  If  there  is  gravity  data  available 
in  the  DOD  Gravity  Library  that  falls  within  the  area,  an  attempt  is 
made  to  evaluate  th^  Ag  values.  This  is  usually  done  by  simple 
comparison  of  the  values.  Many  times,  this  comparison  is  not  possible 
because  of  lack  of  data,  and  the  Ag  is  accepted  simply  because  there 
is  nothing  better. 

6.  Average  of  Smaller  Size  Quads 

In  areas  where  the  gravity  field  has  been  studied  in  detail, 

Agf  values  of  small  size  surface  elements  are  formed.  In  this 
method,  the  Ag^  values  of  the  small  size  surface  elements  are 
averaged  to  form  the  representative  mea:  of  a  larger  surface  element. 
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This  method  Is  always  used  where  smaller  size  means  are  available. 
The  general  equation  Is: 


l 

1=1 


where: 

(igf)1  s  mean  gravity  anomaly  of  the  large  surface  element. 

(£$■ )  =  mean  gravity  anomaly  of  the  small  surface  element. 

n  =  number  of  small  surface  elements  within  the  larje 
surface  element. 

In  Figure  6,  the  block  has  been  divided  ito  5'  x  5'  elements, 
and  the  agf  for  each  5'  x  5'  element  shown.  The  Agf  for  the 
1°  x  1°  1:,  +25  mgals  using: 

144  _ 

(Agf)10  x  s  to  ^Mi 

To  produce  the  basic  5'  x  5'  Agf  values,  an  automated  method 
which  predicts  a  5'  x  5'  Agf  and  Its  accuracy  by  using  point 
gravity  data  Is  used.  The  basic  equation  for  the  least-squares 
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prediction  Is: 


where: 


-1 

C11  C12  *  *  *  Cln 

A  9! 

C21  C22  *  *  *  C2n 

•  • 

A92 

•  %  * 

•  •  * 

• 

• 

Cnl  Cn2  *  *  *  Cnn 

&9n 

Agp  *  predicted  gravity  anomaly  at  a  point 


ag^  =  given  gravity  anomaly 

C] i  *  C22  *  C33  =  variance  of  gravity  anomalies 


=  covariance  of  gravity  anomalies  located  at  points  1  and  k. 


This  equation  Is  derived  and  thoroughly  discussed  by  Moritz  [2]. 
Rapp  [3],  [4],  [5]  developed  a  computer  program  from  these  equations 
which  Is  now  used  to  predict  most  5'  x  5'  mean  anomalies  (either 
Bouguer  or  free-air).  These  references  contain  numerical  examples 
with  discussion  of  the  application  of  this  method. 

There  are  basically  two  methods  of  prediction: 

a.  The  ag^  predicted  directly  from  point  free-air  gravity 
anomalies.  This  method  is  principally  used  to  obtain  predictions 
at  sea  or  In  poorly  mapped  land  areas  where  5'  x  5'  mean  elevations 
are  unavailable. 
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b.  The  Ag^  obtained  from  eleyatton  corrected,  pre^cted 
Agg.  The  Agg  predictions  are  based  on  Bouguer  point  gravi^  anomalies. 
The  S’  x  5'  mean  elevations  are  multiplied  by  the  constant  0.1119 
mgal /meter  and  added  to  the  predicted  Agg  to  obtain  the  predicted 
Ag^.  This  method  is  always  used  in  land  areas  where  5'  x  5'  mean 
elevations  are  available. 

Two  sets  of  covariances  are  computed  for  a  5°  x  5°  area.  One 
set  is  for  east-west  profiles  while  the  other  set  is  for  north-south 
profiles.  These  two  sets  are  combined  into  one  system  of  covariances 
by  fitting  the  data  in  a  least  squares  sense  to  a  tenth  order 
polynomial  in  r,  the  separation  of  the  anomalies,  as  described  by 
Rapp  [4],  [5],  This  may  be  written  as: 

C(r)  =  CQ  +  ^r  +  C2r2  +  .  .  .  C^r10. 

If  a  coefficient  set  has  not  been  computed  for  the  area  to  be 
predicted,  a  computed  set  for  a  similar  area  is  substituted.  The 
choice  of  which  coefficient  set  to  use  is  based  on  a  comparison  of 
the  magnitude  of  the  variance  (Cg)  to  the  magnitude  of  the  variability 
of  the  gravity  anomaly  field.  The  Cg  term  produced  by  a  smooth 
gravity  field  is  numerically  smaller  than  the  Cg  term  produced  by 
a  rough  gravity  anomaly  field. 

The  point  gravity  anomalies  used  for  prediction  are  geographically 
separated  by  a  minimum  of  0.4  minutes  in  latitude  and  longitude. 

This  is  necessary  to  prevent  singular  conditions  occurring  in  the 
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matrix  inversion,  Rapp  [3].  Point  gravity  data  can  be  deleted, 
scaled,  or  translated  thus  allowing  for  the  elimination  of  known 
systematic  errors.  A  predetermined  acceptance  or  rejection  criterion 
expressed  as  a  multiple  of  the  standard  deviation  (e.g.,  2 0)  of  a 
single  observation  computed  from  the  Bouguer  or  free-air  point 
gravity  anomalies  is  used  to  eliminate  extreme  point  gravity 
anomaly  values  from  those  points  used  in  the  prediction.  For 
example,  if  ten  points  are  to  be  used  in  the  prediction  and  one 
additional  point  is  to  be  considered,  the  eleven  closest  points 
are  located,  the  standard  deviation  of  a  single  observation  is 
computed  using  either  the  free-air  or  Bouguer  point  gravity  anomalies, 
and  the  point  which  exceeds  the  rejection  criterion  is  rejected. 

Where  one  or  more  of  the  elevations  of  the  closest  points  is  less 
than  or  equal  to  zero,  the  error  and  rejection  criteria  are  computed 
using  the  free-air  point  gravity  anomalies.  In  all  other  cases,  it 
is  computed  using  the  Bouguer  point  gravity  anomalies. 

After  the  5'  x  5'  Ag^  values  have  been  predicted,  they  are 
placed  in  a  master  file.  The  values  are  than  averaged  to  produce 
the  1°  x  1°  Agf. 

This  method  is  considered  to  be  the  best  for  5'  x  5'  prediction, 
provided  data  to  compute  the  necessary  statistical  coefficients  is 
available.  The  gravity  data  it  areas  where  this  method  is  used 
should  have  a  fair  distribution  and  can  have  any  random  pattern. 

The  method  works  equally  well  in  ocean  and  land  areas. 
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The  1°  x  1°  values  determined  by  this  method  replace  all 
other  Ag.f  values  previously  computed  in  an  area.  The  disadvantage 
in  using  this  method  is  the  slow  process  of  determining  the  Ag^ 
values  for  the  small  size  surface  elements. 

7.  Modified  Average  Free-Air 

This  is  an  automated  method  which  is  used  on  a  worldwide  basis 
to  compute  1°  x  1°  Ag^  values.  In  each  1°  x  1°  area,  the  simple 
average  gravity  anomaly  for  each  10'  x  10‘  area  which  contains 
point  gravity  data  is  computed.  The  10'  x  10'  gravity  anomalies 
are  then  averaged  to  obtain  30'  x  30'  mean  gravity  anomalies. 

These  30'  x  30'  mean  gravity  anomalies  are  averaged  for  the  1°  x  1° 
modified  average  gravity  anomaly. 


(A9f)10'  x  10'  "  n  J  ^A9f^ 

i=l  1 


(A9f*30'  x  30'  ~  k  ^>10'  x  10' 


(A9fV  x  1°  j  J/^f^O*  x  30' 


where: 

n  =  number  of  gravity  anomalies  in  a  10'  x  10'  square, 
k  =  number  of  non-empty  10'  squares  in  a  30'  square, 
j  =  number  of  non-empty  30'  squares  in  a  1°  square. 

In  Figure  7,  the  1°  x  1°  block  has  been  divided  into  10'  x  10' 
blocks,  and  the  point  gravity  anomalies  have  been  plotted  and 
contoured  at  a  10  mgal  interval. 
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Figure  7.  Point  Free-Air  Gravity  Anomalies 


Figure  8  shows  the  simple  average  for  each  10'  x  10'  block  using 
the  data  In  Figure  7. 


Figure  8.  10'  x  10'  Mean  Free-Air  Gravity  Anomalies 
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The  10'  x  10'  blocks  in  Figure  8  are  averaged  for  the  30'  x  3Q' 
blocks  shovn  in  Figure  9. 


Figure  9.  30'  x  30'  Mean  Free-Air  Gravity  Anomalies 
The  modified  average  value  for  the  1°  x  1°  block  is: 


T  |  (4!f)30l  *  30‘  =  *  (89)  =  22  m9al5- 

The  example  can  be  used  to  illustrate  how  the  distribution  of 
the  point  gravity  data  and  the  method  of  computation  can  affect 
the  final  ag^  of  the  block.  The  simple  average  of  the  points  in 
Figure  7  is  12  mgals  which  is  too  small  a  value  for  tne  mean. 

When  using  the  simple  average,  too  much  weight  is  given  to  the 
densely  surveyed  area  in  the  lower  right  corner.  The  average  of 
the  IQ*  x  10'  blocks  in  Figure  8  is  20  mgals.  Here,  the  mean  is 
too  low  due  to  the  distribution  of  the  10'  x  10'  blocks  which 
contain  gravity  data.  By  comparing  the  computed  modified  average 
(22  mgals)  with  the  contours  of  the  example  (Figure  7),  the  modified 
average  best  describes  the  ag^  of  the  block. 
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Satisfactory  results  are  usually  obtained  from  the  modified 
average  method.  In  areas  where  the  distribution  is  very  poor 
and/or  where  the  point  count  Is  less  than  5,  the  computed  1°  x  1° 

Ag.p  may  have  to  be  adjusted  In  order  to  be  more  compatible  with 
the  physical  situation. 

No  Ag^  values  are  extrapolated  or  interpolated  for  the  10'  x  10' 
blocks  which  do  not  contain  gravity  data.  Data  which  is  immediately 
adjacent  to  the  1°  x  1°  block  is  not  used,  and  no  mathematical 
correlation  to  the  1°  x  1°  blocks  surrounding  the  area  is  attempted 
here.  The  1°  x  1°  Ag^  values  are  visually  inspected  to  insure  that 
they  are  generally  compatible  with  surrounding  values  and  the  physical 
conditions  of  the  area. 

This  method  is  used  to  rapidly  incorporate  point  gravity  data 
Into  the  worldwide  1°  x  1°  gravity  anomaly  file.  The  modified 
average  values  are  replaced  in  areas  where  the  gravity  field  is 
studied  more  intensely  using  other  computational  techniques. 

Most  1°  x  1°  Agf  values  have  been  computed  using  the  Modified 
Average  Free-Air  method.  The  approximate  number  and  percentage 
of  the  1°  x  1°  Agf  values  computed  by  the  various  methods  are 
tabulated  in  Table  1. 
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TABLE  1 


NUMBER  AND  PERCENTAGE  OF  WORLDWIDE  1°  x  1°  Agf  VALUES 
COMPUTED  BY  EACH  METHOD  T 


METHOD 

NUMBER 

PERCENTAGE 

Bouguer  Gravity  Anomaly  Map  Estimates 

857 

3.1 

Free-Air  Anomaly  Map  Estimates 

692 

2.5 

Geophysical  Correlation 

5691 

20.7 

Interpolation  or  Extrapolation 

228 

1.0 

Mean  Anomalies  Received  from  Other  Organizations 

2045 

7.4 

Average  of  Smaller  Size  Quads 

948 

3.4 

Modified  Average  Free-Air 

16980 

61.9 

TOTAL 

27441 

100.0 
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COMPUTATIONAL  METHODS  FOR  DETERMINING 


THE  ACCURACY  OF  A  MEAN  FREE-AIR  GRAVITY  ANOMALY 


GENERAL  ERROR  EQUATION 


The  standard  error  (a _ }  of  a  computed  or  predicted  mean 

Ag 

gravity  anomaly  is  in  general  determined  by: 


where: 

Ep  =  error  in  the  point  gravity  anomalies  which  are  used  in 
the  mean  gravity  anomaly  (Ag")  computation  or  prediction. 

Er  =  error  in  the  mean  elevation  value  (land), 
h 

Er  *  error  of  representation  which  is  a  function  of  the  density 
and  distribution  of  th_  point  anomalies  and  the  gravity 
gradient. 

1.  Error  in  the  Point  Gravity  Anomalies 

The  master  point  gravity  anomaly  file  contains  both  evaluated 
and  unevaluated  gravity  data.  Surveying  organizations  transmit 
gravity  data  which  may  or  may  not  contain  an  estimate  of  the  survey 
error  (Ep).  When  gravity  data  is  received,  it  usually  is  immediately 
placed  in  the  &g  file.  After  being  placed  in  the  Ag  file,  the  source 
may  be  evaluated  and  assigned  an  Ep  value,  or  a  previously  determined 
Ep  value  may  be  changed. 

When  computing  or  predicting  Ag^  values  using  an  evaluated  source, 
Ep  is  usually  obtained  directly  from  the  Ag  file.  When  using  an 
unevaluated  source,  Ep  is  determined  empirically  by  comparing  Ag 
values  of  the  unevaluated  source  with  Ag  values  of  adjacent  source(s). 
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Usually*  a  computer  print  plot  of  the  Agg  values  is  used  for 
comparison  purposes.  Since  the  Ag  file  contains  both  evaluated 
and  unevaluated  sources,  the  empirically  determined  Ep  is  used 
when  an  unevaluated  source  is  encountered.  If  sufficient  evaluated 
data  is  available  for  prediction,  the  unevaluated  data  is  deleted. 

In  computing  or  predicting  Ag^  values,  Ep  is  always  set  equal  to 
or  greater  than  one  milligal  every^hough  the  source  is  evaluated 
to  be  better  than  one  milligal.  This  is  because  the  smallest  final 
accuracy  recorded  for  a  Ag^  value  is  +1  mgal . 

2.  Error  in  the  Mean  Elevations 

Mean  elevation  errors  (E-)  are  obtained  directly  from  the 

master  mean  elevation  file.  The  mean  elevation  accuracies  are 

determined  at  the  time  the  mean  elevations  are  determined.  The 

mean  elevation  error  is  converted  to  mill i gals  (Er  =  0. » 11 9  x  E-). 

h  e' 

When  the  Agf  is  computed/predicted  directly  on  the  geopotential 
surface  (ocean  area),  Ag^  values  are  used  directly.  Normally,  the 
stated  accuracies  of  the  mean  elevation  values  are  accepted.  However, 
if  there  are  indications  that  a  mean  elevation  or  its  stated  accuracy 
is  questionable,  the  mean  elevation  or  its  stated  accuracy  is  changed 
accordingly. 

3.  Error  Representation  (Prediction) 

Er  is  the  error  of  a  computed  or  predicted  mean  gravity  anomaly 
without  any  consideration  for  errors  due  to  inaccuracies  in 
Ag  values,  mean  elevations,  or  similar  quantities.  It  is 
dependent  only  upon  the  density  and  distribution  of  the  point 
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gravity  anomaly  data  used  to  compute  or  predict  the  mean  gravity 
anomaly  and  the  gradient  of  the  gravity  anomaly  field. 

ACCURACY  COMPUTATION 
1 .  Bouguer  Gravity  Anomaly  Map  Estimates 

The  accuracy  determination  for  this  method  is  subject  to  many 
variables  which  in  most  cases  can  only  be  estimated.  These 
variables  include: 

Accuracy  of  Bouguer  Gravity  Anomaly  Map 
Accuracy  of  point  source 
Accuracy  of  map  compilation 
Accuracy  of  Estimation 

Distribution  and  density  of  point  control 
Gravity  anomaly  field  gradient 
Quad  size 

The  accuracy  of  the  map  depends  on  the  accuracy  of  the  point 
gravity  data  used  for  control,  and  the  method  and  care  used  in 
producing  the  map.  In  most  cases,  these  accuracy  values  are  not 
indicated  on  the  maps  and  are  not  obtainable  from  other  sources. 
Therefore,  the  analyst  must  estimate  some  accuracy  value  for  the 
map  based  mainly  on  his  experience.  Evaluated  data  can  sometimes 
be  used  to  evaluate  a  map  if  the  data  lies  within  or  adjacent  to 
the  area. 

For  the  estimation  of  the  accuracy  of  the  estimated  Fg  ,  the 

B 

analyst  has  a  little  more  control.  If  the  points  are  plotted, 
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their  distribution  and  density  will  help  In  the  accuracy 
estimation.  If  they  are  not,  the  detail  in  the  contours  may 
indicate  point  distribution.  The  contribution  to  the  error  of 
estimation  by  the  gradient  of  the  field  must  be  considered  in 
conjunction  with  the  quad  size.  This  error  of  estimation  is 
made  at  the  same  time  as  the  mean  gravity  anomaly  is  estimated, 
and  there  is  no  general  equation  used.  The  root  mean  square  of 
the  Bouguer  gravity  anomaly  map  accuracy  and  the  accuracy  of 
estimation  is  the  final  accuracy  for  the  1°  x  1°  quad. 

2.  Free-Air  Gravity  Anomaly  Map  Estimates 

The  accuracy  determination  for  this  method,  like  that  for 
the  Bouguer  Map  Estimation  Method,  is  subject  to  many  variables 
which  in  most  cases  can  only  be  estimated. 

The  general  problems  here  are  the  same  as  those  in  the 
Bouguer  Map  Estimation  Method  discussed  in  Part  II,  Section  1. 

The  major  difference  here  is  the  much  steeper  gradients,  and 
the  inherent  difficulty  in  estimating  the  contribution  of  this 
steeper  gradient  to  the  accuracy  estimation. 

3 •  Geophysical  Correlation 

The  many  different  methods  used  to  compute  the  accuracies  of 
geophysically  correlated  values  are  recorded  in  the  document 
referenced  in  Part  I,  Section  3. 

4 .  Interpolation  or  Extrapolation 

The  accuracy  of  an  interpolated  or  extrapolated  gravity 
anomaly  is  at  best  an  estimation  made  by  the  analyst  based  on  the 
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particular  situation  and  his  experience.  The  amount,  distribution, 
and  quality  of  observed  gravity  data,  and  the  bathymetric  or 
topographic  characteristics  of  the  area  are  considered.  In  many 
cases,  a  blanket  error  is  assigned  to  areas  of  rather  large  extent. 

5.  Mean  Anomalies  Received  from  Other  Organizations 

The  methods  and  data  used  by  other  organizations  to  compute 
Ag  values  and  their  related  accuracies  often  are  not  documented  or 
available.  If  evaluated  gravity  data  in  the  DOD  Gravity  Library  lies 
within  the  area,  an  attempt  is  made  to  validate  the  accuracy  values 
furnished  by  the  source.  Usually  this  is  not  possible  because  of  the 
lack  of  gravity  data,  and  the  furnished  accuracy  values  are  recorded. 

6.  Average  of  Smaller  Size  Quads 

The  general  form  of  the  empirical  equation  for  the  accuracy  of 
a  Agf  made  up  of  the  average  of  smaller  quads  is: 

r  p-,1/2 

KinT 

where: 

a.j  =  error  of  large  surface  element. 
o.j  =  error  of  small  surface  element. 

N  =  number  of  small  surface  elements, 
n  *  number  of  small  surface  elements  with  gravity  data. 

In  Figure  10,  the  1°  x  1°  block  has  been  divided  into 
5'  x  5'  elements,  and  the  accuracy  of  the  mean  gravity  anomaly 
for  each  5*  x  5'  element  is  shown. 
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Figure  10.  Accuracy  of  5'  x  5'  Mean  Gravity  Anomalies 
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The  a-j  for  the  1°  x  1°  is  +1.7  mgals  using: 


This  accuracy  would  then  be  recorded  as  +2  mgals.  Since  the 
5‘  x  5'  Ag^.  values  have  been  studied  in  detail  and  edited,  the 
1°  x  1°  Ag^.  values  and  their  accuracies  can  be  considered  the 
best  available. 

The  basic  equation  for  the  mean  square  error  of  least  squares 

prediction  is  used  to  compute  the  error  of  the  5'  x  5'  Ag  (Moritz  [2]). 

f 


where: 

Cq  a  C  =  C  =  variance  of  gravity  anomalies 

=  covariance  of  gravity  anomalies  located  at  points  i  and  k. 
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The  localized  error  covariances  are  discussed  In  Part  I, 

Section  6.  The  computation  of  the  covariances  Is  fairly  straight¬ 
forward,  though  lengthy.  The  specific  computations  needed  for  the 
covariance  computation  have  been  fully  described  In  Rapp  [3],  [4], 
and  [5].  The  screening  and  elimination  of  erroneous  point  gravity 
data  and  methods  to  eliminate  systematic  errors  are  addressed  in 
Part  I,  Section  6.  This  method  is  considered  the  best  and  is  used 
to  predict  most  5'  x  5'  values  and  their  accuracies. 

7.  Modified  Average  Free-Air 

In  the  computation  of  the  standard  deviation  of  the  modified 
average,  10'  x  10'  blocks  with  gravity  data  are  considered  independent, 
and  those  10'  x  10'  blocks  without  gravity  data  are  considered 
dependent.  The  general  form  of  the  equation  is: 


'1°  x 


r  -  1  (n(RMS)J  2  +  (36-n)  (M)2 
36  l—1  J 


1/2 


where: 

RMS  = 


n 

1 

i=l 

m2 

L  n 

1/2 


m  =  a  dependent  error  =  error  of  a  blank  10*  x  10'  block  = 
+19  mgal. 

M  =  an  independent  error  =  error  of  a  10'  x  10'  block  with 
gravity  data  =  +2  mgal. 
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n  s  number  of  10'  x  10'  blocks  without  gravity  data. 

Using  these  values,  the  preceding  general  equation  reduces  to: 

r  1  ^2 

c  •  1  (n-19)2  +  (36-n)  (2)2 

1  x  l  36  L  J 

In  Figure  11,  the  10'  x  10'  blocks  marked  with  X  contain  point 
gravity  data.  The  number  of  blocks  without  point  data  is  25. 


Figure  11.  10'  x  10'  Blocks  Containing  Point  Gravity  Data 

1/2 

0  =  f[25(19)]2  +  (11)(2)2  ] 

10  *1°  36 


1°  x  1 


[  225625  +  44  1 
36 


1/2 


V  x  1° 


+13  mgals. 
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Since  there  are  no  Ag^  values  extrapolated  or  interpolated 
for  the  blank  10'  x  10'  blocks,  the  errors  for  these  blank  blocks 
are  likewise  not  extrapolated  or  interpolated.  Instead,  a  blanket 
value  of  +19  mgals  is  used.  This  value  was  taken  from  Heiskanen 
and  Moritz  [1],  pg  279,  who  stated  that  the  error  variance  for  a 

o 

1°  x  1°  block  with  one  gravity  station  is  356  mgals  .  The 
independent  blocks  have  been  assigned  a  blanket  standard  deviation 
of  +2  mgals.  Again,  this  method  was  developed  in  order  to  produce 
a  1°  x  1°  working  field  from  worldwide  point  data.  As  areas  are 
more  intensely  studied,  these  modified  averages  are  replaced. 

SUMMARY 

This  publication  describes  the  methods  presently  used  by  DMAAC 
to  determine  mean  gravity  anomalies  and  their  accuracies.  The 
definition,  equations,  sample  computation,  advantages,  and  disadvan 
tages  of  each  method  are  discussed. 

Research  to  find  better  methods  to  compute  mean  anomalies  and 
their  accuracies  is  a  continuing  assignment.  Some  of  the  equations 
now  in  use  have  been  empirically  derived.  As  research  progresses, 
the  current  methods  may  be  discontinued  or  altered  to  incorporate 
the  newly  developed  improvements. 
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